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Abstract. In this paper we develop a quantitative version of Enss' method to 
establish global-in-time decay estimates for solutions to Schrodinger equations 
on manifolds. To simplify the exposition we shall only consider Hamiltonians of 
the form H := — ^ Aj\,/, where is the Laplace-Beltrami operator on a mani- 
fold M which is a smooth compact perturbation of three-dimensional Euclidean 
space which obeys the non-trapping condition. We establish a global-in- 
time local smoothing estimate for the Schrodinger equation ut = —iHu. The 
main novelty here is the global-in-time aspect of the estimates, which forces 
a more detailed analysis on the low and medium frequencies of the evolution 
than in the local-in-time theory. In particular, to handle the medium frequen- 
cies we require the RAGE theorem (which reflects the fact that H has no 
embedded eigenvalues), together with a quantitative version of Enss' method 
decomposing the solution asymptotically into incoming and outgoing compo- 
nents, while to handle the low frequencies we need a Poincare-type inequality 
(which reflects the fact that H has no eigenfunctions or resonances at zero). 



1. Introduction 



Let {M,g) = (R^,!?) be a compact perturbation of Euclidean space^ R'^, thus M is 
R^ endowed with a smooth metric g which equals the Euclidean metric outside of 
a Euclidean ball B{0,Ro) {x e R^ : < i?o} for some fixed Rq- We consider 
smooth solutions to the Schrodinger equation 

(1) Ut = —iHu 

where for each time t, u{t) : M — > C is a Schwartz function on M, and H is the 
Hamiltonian operator 

where A^/ ^"'Vj is the Laplace-Beltrami operator (with denoting covariant 
differentiation with respect to the Levi-Civita connection, in contrast with the 
Euclidean partial derivatives dj). Note that H is positive definite and self-adjoint 
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^The analysis we give here also extends to higher dimensions n > 3, and are in fact slightly 
easier in those cases, but for simplicity of exposition we restrict our attention to the physically 
important three-dimensional case and to compact perturbations. The hypothesis that M is topo- 
logically is technical but in any event is forced upon us by the non-trapping hypothesis which 
we introduce below; see 1141 . 
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with respect to the natural inner product 

(2) {u,v)l2^m) ■= / u{x)v{x) dg{x), 

Jm 

where dg := -^/det gij(x)dx is the standard volume element induced by the metric 
g. In fact the spectrum of H consists entirely of absolutely continuous spectrum on 
the positive real axis [0, +00), in particular H has no eigenvalues or resonances at 
any energy. In particular H enjoys a standard functional calculus on L^{AI), the 
Hilbert space associated to the inner product, and one can define the homogeneous 
Sobolev norms ||u||^s(jv/) '■~ II-^*^^"II-L^(a/) on Af using fractional powers of H for 
all — 1 < s < 1. (One can of course define these norms for other s also, but there 
are some technicalities when s is too negative that we do not wish to address here) . 
Thus for instance ||w||^o(^,/) = II'"IIl2(m) and 

(3) \\u\\\,^j^_j^ = {u,Hu)l2(m) dg. 



It is well known (see [2]) that for any time to and any Schwartz initial data uq, 
there exists a unique global-in-time Schwartz solution u : R x ill — > C to 
with initial data 11(0) = uq, indeed we have u{t) — e^**^uo. In this paper we 
develop a quantitative variant of Enss' method to obtain a new global-in-time local 
smoothing estimate for such solutions; to avoid needless technicalities we shall 
always restrict ourselves to Schwartz solutions. The methods here extend to more 
general classes of Hamiltonians than those considered here, for instance they can 
handle asymptotically flat manifolds of dimension n > 3 as well as short-range 
potentials provided that there are no resonances or eigenfunctions at zero; however 
to simplify the exposition we have chosen to restrict attention to the simple case of 
zero potential and compact perturbations of three-dimensional Euclidean space. We 
shall pursue more general Hamiltonians in |5] using a somewhat different method 
(based on limiting absorption principles). 

Let us begin by recalling some earlier results. In the case where M is Euclidean 
space R", so that H = Hq :— — ^Ap^a is just the free Hamiltonian, then we have 
the well-known global-in-time local smoothing estimate (see QI], which 
we shall phrase as 
(4) 

/OO 

for any a > and any Schwartz initial data /, where (x) := (1 + jxp)^/^. It is well 
known that the condition on a is sharp (there is a logarithmic divergence in the 
left-hand side when cr = 0). 

Now we suppose that {M,g) = (R^,^) is a compact perturbation of Euclidean 
space obeying the non-trapping condition. It is known (see 0) that one has 
the local-in-time local smoothing estimate 
(5) 

||(x)-l/2-'^Ve-^*^^.o|li.(M) + ll(^)-'/'-'^e-^*^«o|li.(M) dt < C,.IMho\\l^,.^M) 
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for all compact time intervals / C R, where the constant on the right-hand side 
is allowed to depend on /, if and only if the manifold M is non-trapping, i.e. 
every geodesic s <—>■ x{s) in M eventually goes to spatial infinity as have a theory 
for the Schrodinger equation as s ^ ±cx). To understand why the non-trapping 
condition is necessary, observe that the localization in time means that the low and 
medium energies are easily controlled (for instance by using the Sobolev embedding 
H^^^{AI) C L'^[M) and Holder's inequality, noting that at low and medium energies 
yg-jtHy^ is also in H^^^{AI)), and so one only needs to understand the evolution 
of the high energies, which evolve semi-classically. The semi-classical limit of the 
estimate © is the estimate Jji{x{s))~^~'^'^ \i{s)\'^ ds < Ccr,M\i{0)\ for any geodesic 
s x(s), which is easily seen to hold if and only if the manifold is non-trapping 
(this also explains the requirement that ct > 0). Of course, one needs semi-classical 
tools such as pseudo-differential operators and the positive commutator method in 
order to make this argument rigorous; see j2|, |S1 (or SectionlHl for more details. 

The main result of this paper is to unify the global-in-time Euclidean estimate (0J 
with the local- in-time manifold estimate © as follows. 

Theorem 1.1. Let M be a smooth compact perturbation of R"^ which is non- 
trapping and which is smoothly dijfeomorphic to R'^ . Then for any Schwartz solu- 
tion uit, x) to and any cr > we have 
(6) 

/ ||(x)-V2-.Ve-*^uo||i.(M) + ||(x)-3/^-'^e-*^«o|li.(M) dt < C^M^oWl^;.,^)- 
JR 

In other words, the constant Ca-.i,M in ® can be taken to be independent of the 
interval I . 



As mentioned earlier, this is not the strongest result that one could obtain with 
this method; our purpose here is merely to illustrate a model example in which the 
method applies. By using different methods (in particular the limiting absorption 
principle) we were able to obtain results for more general Hamiltonians, see [S]. 
However we believe the method we present here still has some merit, in particular 
it is conceptually straightforward and seems to have some hope of generalizing to 
more "time-dependent" or "non-linear" situations in which spectral theory tools are 
less useful. For instance, the methods here were already used in to obtain new 
results on the asymptotic behavior of focusing non-linear Schrodinger equations, at 
least in the spherically symmetric case. 

We now informally discuss the proof of the theorem. First note that this estimate is 
already proven when M is Euclidean space R"^, and we are only considering mani- 
folds M which are compact perturbations of Euclidean space. Thus it is reasonable 
to expect that the only difficulties in proving (j^ will arise from the compact region 
B{0, Rq), and indeed by modifying the proof of I^J we will be able to show that a 
global-in-time local smoothing estimate in a slight enlargement -6(0, R) of -6(0, Rq) 
will automatically imply the full estimate ((HJ- Thus we may (heuristically, at least) 
restrict our attention to a ball such as .8(0, R). 
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Next, we use the spectral theorem to decompose the evolution into low-energy, 
medium-energy, and high-energy components. The high energy components turn 
out to be treatable by the same positive commutator arguments used to prove 
the local-in-time estimate (jSj). From a heuristic viewpoint, this is because high- 
energy components propagate very quickly, and thus only linger in the compact 
region i?(0, R) for a very short period of time, after which they escape into the 
Euclidean region of M and never return to -6(0, i?) again. Of course the non- 
trapping hypothesis is essential here. At a more technical level, the reason why we 
can adapt the arguments used to prove © is that the error terms generated by the 
positive commutator method are lower order than the main term, and can thus be 
absorbed by the main term in the high-energy regime even when the time interval 
/ is unbounded. 

The low-energy components are easy to treat, but for a different reason, namely that 
there is an "uncertainty principle" that shows that solutions which have extremely 
low energy cannot be concentrated entirely in the compact region K, and thus the 
low energies cannot be the dominant component to the local smoothing estimate 
inside this region K. This is ultimately a reflection of the well-known fact that a 
Hamiltonian H = —^Ai\i with no potential does not have any resonances or bound 
states at zero or negative energies. 

The most interesting component to treat is the medium energies regime, which 
requires new methods. These are energies which are not high enough to behave 
semi-classically and escape the compact region K by means of the non-trapping 
hypothesis, but which are not low enough to escape the compact region K by 
means of the uncertainty principle. Fortunately there is a third mechanism by 
which we can force solutions of the Schrodinger equation to escape the compact 
region K, namely the Ruelle-Amrein-Gorgescu-Enss (RAGE) theorem. The point 
is that (as is well known) H contains no embedded eigenfunctions in the medium 
energy portion of the spectrum (or indeed anywhere in the spectrum), and so the 
RAGE theorem then ensures that any given solution to |^ must eventually vacate 
the region K after some time T. We recall the abstract version of the RAGE 
Theorem (see e.g. jH]). 

Theorem 1.2 (RAGE). Let H be a self-adjoint operator on a Hilbert space TL. If 
C :T-C ^ H is a compact operator and uq lies in the continuous subspace of H , then 



The RAGE theorem played a crucial role in Enss' approach to scattering for 
Schrodinger operators H = Hq + V, In fact our treatment of medium en- 
ergies can be viewed as a quantitative version of the Enss' method in a sense that 
we use both the RAGE type inequality and the decomposition into incoming and 
outgoing waves, also a major part of the Enss' work, but derive an a priori estimate, 
as opposed to a qualitative result about completeness of the wave operators. The 
other major difference is that in scattering theory, density arguments allow one to 
consider only medium energies and compactly supported data, whereas we must 
necessarily treat all energy ranges and allow our data to have arbitrary support. 
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One might object that the RAGE theorem is "quahtative" in nature, in that the 
time T required for a solution to leave K depends on the choice of solution and 
thus need not be uniform. However because we have localized the solution in both 
frequency (to medium energies) and position (to the region K), and because we can 
use linearity to normalize the H^^^{M) norm of u, the solution is in fact effectively 
contained in a compact region of (quantum) phase space. Because of this, one can 
make the time T required for a solution to leave K to be uniform for all medium- 
energy solutions u. 

There is however still a remaining difficulty for medium energy solutions, which is 
that once a solution leaves the compact region K one needs to ensure that it does 
not return back to K, since a solution which periodically left and then returned to 
K would eventually contribute an infinite left-hand side to To resolve this we 
introduce a quantitative version of Enss' method. The starting point is the obser- 
vation that, in the exterior of the domain K, any function can be decomposed into 
"outgoing" and "incoming" components (very roughly speaking, this corresponds to 
the spectral projections X[o.+oo){'i9r) and X{-oo,o]{'^9r) where dr is the radial deriv- 
ative, although we shall not perform these projections directly due to the singular 
behavior of the operator idr). Outgoing components will evolve towards spatial 
infinity as t +00, whereas incoming components will evolve towards spatial in- 
finity as t —^ —00. In particular in both cases the solution will not encounter the 
compact region K and the evolution is essentially Euclidean in nature. 

Now suppose that a medium energy solution to (QJ is localized to K at some time 
to- By the quantitative RAGE theorem, at some later time to + T the solution has 
mostly vacated the region K. By Enss' decomposition it can have either outgoing 
or incoming components. But one can show that there is almost no incoming 
component, because if we evolved backwards in time from tg + T back to to we see 
that the incoming component would have evolved back to a region far away from K , 
and thus orthogonal to the initial data. Thus at time to + T, the solution consists 
almost primarily of outgoing components. But then by the preceding discussion this 
means that the solution will continue to radiate to spatial infinity for times after to, 
and thus never return to K. To summarize, we have shown that any component of a 
medium energy solution which is located in K at time to, will eventually radiate to 
spatial infinity as i — > -l-oo; a similar argument also handles the t — > —00 evolution. 
Combining this with the finite energy of u (note that in the medium energy regime 
all Sobolev norms are equivalent), one can obtain the estimate 

The global- in-time local smoothing estimate in Theorem II . 1 1 has a number of con- 
sequences; for instance by combining it with the arguments of Staffilani and Tataru 
|12j one can obtain global-in-time Strichartz estimates for compact non-trapping 
perturbations of Euclidean space, which then can be used to transfer some local and 
global existence results for non-linear Schrodinger equations on Euclidean space, to 
the setting of compact non-trapping perturbations. We will not discuss these (fairly 
standard) generalizations here, but see ^ for further discussion. 
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Throughout this paper, the manifold M , the radius i?o, and the exponent a will be 
fixed. All constants C will be allowed to depend on a and M (and hence on Rq). 
If C needs to depend on other parameters we will indicate this by subscripts. 

Suppose Aq, Ai, . . . , Afe are real parameters. We use OAi_,...,Ak{^) to denote any 
quantity depending on Ai,...,Ak and possibly some other quantities, which is 
bounded in magnitude by some constant CAi,...,Ak- For some fixed c (usually c — 
or c = oo), we also use OAo^c;Ai,...,Aki^) to denote a quantity depending on 
Aq, Ai, . . . , Ak and possibly some other parameters, which is bounded in magnitude 
by some quantity FAi,....A^iAo) such that 



Often k will be zero, in which case the above notations simply read 0(1) and 
Oyio— c(l) respectively. We also abbreviate OAi,...,Ak{^)^ and OAo^c;Ai,...,Aki^)^ 
as OAi,...,AkiX) and oao^c;Ai....,a^{X) respectively. 

In the absence of parentheses, we read operators from right to left. Thus for in- 
stance, Vfg denotes the function V{fg) rather than {V f)g. 

We use the usual summation conventions on indices, and use g^^ to denote the dual 
metric to gjk on the cotangent bundle, and use g to raise and lower indices in the 
usual manner. We use Vj, V-' to denote the usual covariant derivatives with respect 
to the Levi-Civita connection on M; these can be applied to any tensor field; we 



use |V/|g — y'^^f'^jf to denote the magnitude of the gradient with respect to 
the metric g, and |V/| to denote the Euclidean magnitude of the gradient. Since 
Vg — 0, the indices of covariant derivatives can be raised and lowered freely, thus for 
instance Am = V-' Vj- = VjV-' . Also these covariant derivatives are anti-selfadjoint 
with respect to the inner product Q, and thus we can integrate by parts using 
these derivatives freely. 

As the operator H is self-adjoint and has spectrum on [0, cxd) we can construct spec- 
tral multipliers f{H) for any measurable function / : [0, oo) ^ C of at most poly- 
nomial growth; in particular we can define fractional powers 

as well as Schrodinger propagators e~**-^ and Littlewood-Paley type operators on 
H. These spectral multipliers commute with each other, and are bounded on if 
their symbol / is bounded. 



2. Notation 



Ao^c 



lim -F!Ai,...,Afc = for all choices of ^i, . . . , Ak- 
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3. Overview of proof 



We now begin the proof of Theorem ll.il The first step wiU be to show that one 
can freely pass back and forth between the slowly decaying weight (a;)~^/^~°' in © 
and a suitably chosen compactly supported weight Lp. By shrinking a as necessary 
we may assume that < cr ^ 1. 

Fix any compact time interval [0,T], and let K{T) be the best constant for which 
the inequality 
(7) 



holds for any Schwartz function uq. From the local-in-time local smoothing theory 
in j2j, |5| we already know that K{T) is finite for each T. Our task is to show that 
K{T) is bounded independently of T; the negative times can then be handled by 
time reversal symmetry. 

Recall that M is equal to Euclidean space in the exterior region \x\ > Rq. To take 
advantage of this, let us fix : M — > R to be a smooth function which equals 1 
when |x| < 4i?o and equals when |x| > 8Rq. We then define the localized quantity 
K^{T) to be the best constant such that 

(8) £ ll^|Ve-"«7.olslli.(Af) dt < K^iTnuoWl,,,^^,^^ 

holds for any Schwartz solution u to (QJ. It is clear that Kip{T) < CipK{T). In 
Section 01 we shall establish the converse inequality 

(9) K{T)<C^ + C^K^{T) 

for all times T > 0, where the constants depend on Lp but not on T. 

In light of we see that to bound K{T) it suffices to bound the localized quantity 
K^{T). We make the technical remark that K^{T) is only required to control first 
derivatives of e^**^Mo, and not e~'*^wo directly. This will be important in the low 
freqency analysis later on. 

The next step is energy decomposition into the very low energy, medium energy, and 
very high energy portions of the evolution. Let < £o ^ 1 be a small parameter 
(depending on ip) to be chosen later, and decompose 1 — Pio + Pmed + Phi, where 
Pio, Pmed, Phi are the spectral multipliers 

Pio:=xiH/so); Prned^xi£oH)-xiH/eo); P,„ := 1 - x(eoi?) 

and X ■ ^ R- is a bump function supported on [—1,1] which equals 1 on 
[—1/2, 1/2]. We shall prove the following three propositions, in Sections IZl HI 
respectively: 



Proposition 3.1 (Low energy estimate). For any T > 0, we have 

yVPioe-'^'^UoWh^M) dt < Oeo^O;^(if(T)2)||,io|llfV2(M)- 
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Proposition 3.2 (Medium energy estimate). For any T > 0, any time-step t > 1 
and radius R > 10i?o, we have 
T 

Proposition 3.3 (High energy estimate). For any T > 0, we have 

T 

||(^VP^,e-^*^Uo||i.(M) dt < o,„^o-AK{Tf )\\uo\\l,,,,^,y 
^ ' 

Combining these three propositions using the triangle inequality and using the 
definition Q of Kip{T), we see that 

{k{tY) 

for any £o > 0, r > 1, and R>1. If we choose r sufhciently large depending on (p, 
£o and i?, and R sufficiently large depending on Lp^ Sq, we conclude that 

combining this with © we obtain 

letting Eq be sufficiently small depending on Lp and recalling that K{T) is finite, we 
conclude that K{T) < dp for all time T, which gives Theorem ll.il 

It remains to prove © and Propositions 13.11 13.21 13.31 This will be done in the 
following sections. 



4. Physical space localization 

We first prove ©. Let u — e^'*^wo be a Schwartz solution to Let T > 0. 
We will allow all constants C to depend on ip and will no longer mention this 
dependence explicitly. We normalize ||uo|liji/2(j(,f') = 1, which by unitarity of e" 
and spectral calculus implies that 



-itH 



(10) sup ||w(t)||^i/2(M) = 1- 

teR 



Our task is to show that 







||(a:)-i/2--vw||i.(M) + dt < (c + cx^(r)2)||woll^v.(M)- 



Note that it will not be relevant whether we measure the magnitude of Vw using 
the metric g or the Euclidean metric as they only differ by at most a constant. 



From we already have 

(11) / / |Vup dxdt < CK^{Tf 

Jo J\x\<3Ro 
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From an easy Poincare inequality argument we also can show that 

/ \u\^dx<C{( \Vu\^dx+[ |Mp dx. 

J\x\<3Ro J\x\<4Ra J3Ro<\x\<4Ra 

Thus it will suffice to work in the Euclidean region > 3i?o and prove that 

(12) r [ Ixl^^-^^lVul^ + lxl'^-^^lu]^ dxdt<C + CK^{T)^. 

Jq J\x\>3Ro 

We now invoke the positive commutator method. Let A be an arbitrary linear 
operator on Schwartz functions. From the self-adjoint nature of we observe the 

Heisenberg identity 

(13) ^(Aw(t),u(t))i2(M) = {i[H,A]u{t),u{t))L^M) 

where i[H, A] = i{HA — AH) is the Lie bracket of H and A; integrating this in t 
and using finil and the duality of H^/^{M) and H-'^/^{M), we obtain 

(14) I ^ {t[H, A]u{t), u{t))mM) dx\ < C\\A\\ 

The positive commutator method is based on choosing A so that so that i[H,A] is 
mostly positive definite, in order to extract useful information out of H14|l . 

Let us now set A equal to the self-adjoint first-order operator 

A ~ —ia ^dk — idkOi,k 

where a : -Af ^ R is the function a xd^^l ~ where x is a smooth cutoff 

supported on the region > 2i?o which equals 1 when |a;| > 3i?0: and < £ ^ 1 is 
a sufficiently small constant depending on i?Oi a,fc denotes the Euclidean derivative 
of a in the direction, and we are summing indices in the usual manner. Since 
Va = 0(1) and V^a = 0(1/|2;|), we observe from Q and the classical Hardy 
inequality Hm/I^^I lli2(R,-'') < ^ll"lli?i(R^) ^^^^ ^ maps H^{M) to L'^{M), and by 
self-adjointness also maps L'^{M) to H~^{M). By interpolation we conclude that 
A maps ij^/^ to H^^^^. Also, since H = —^djdj on the support of x, we can 
compute 

i[H, A] = -^[djdj,a^kdk + dka,k] 

= ^'^{'^Mjj^k + '2a,kjdjk + dkO^kjj + '^dka^kjdj) 

= —2dja^kjdk — 2'^,jjkk 
and thus from (|14|l and an integration by parts we conclude that 

(15) \[ [ 2ajkdkudfu--A'L3a\u\'^dxdt\<C. 

Jo J\x\>Ro ' 2 ti. 

Let us first consider the portion of 115() on the region |a;| > 3Rq, for which a — 
|a;| — Then a computation shows that a has some convexity if e is sufficiently 
small, indeed in this region we have 

|V-u|2 

ajkdkudjU > Ce ; A^a < -cj\x\^+'' 
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for some > 0. Invoking (|15|l and using to estimate the region where 2Ro < 
\x\ < 3Rq, we conclude that 



2 r / riT^+Cer^ dxdt < C,+C,K^iTf+^ f [ Ai^salu]^ dxdt. 

To conclude ((T^ . it thus suffices by ((TT)l to establish following the fixed-time esti- 
mate: 

Lemma 4.1. If e > is sufficiently small, then there exists a constant C > such 
that 



f A|^3a|/P<c/ 

J\x\<3B.o J\x 

for all smooth functions f. (Note that the left-hand side can be negative). 



\x\<3B.o 



Proof. It is possible to establish this from Poincare inequality and the Green's 
function computation below (the main point being that A^a has negative mean), 
but we shall use a compactness argument instead. Let i5 > be a small number to 
be chosen later. It clearly suffices to show that there exists C > such that 



{6 + Al,sa)\f\^<S \f? + C |V/p 



T?' 

\x\<3Ra J\x\<3Ro J \x\<3Ra 

for all smooth functions /. We have chosen to use the Euclidean measure dx here 
but one could equally well run the following argument using the measure dg. 

Suppose for contradiction that the above estimate failed. Then we can find a 
sequence /„ of smooth functions with the normalization 



/ \fn\'+n[ \Vfn\' = l 

J\x\<3Ro J\x\<3Ro 



such that 



limsup / {6 + A^a)\fn\^ > 1. 

n^oo J\x\<3Rn 



n^oc J|a:|<3flo 

By Rellich compactness we can find a subsequence fn- of /„ which converges in 
to a limiting object / G L^{B{0, 3Ro)), and then by Fatou's lemma 

/ |v/p = o. 

J\x\<3Ro 

In other words, / is equal to a constant on the region (x) < SRq. But by Green's 
formula we have 

/ S+A'^a = OiS)+ [ -^Aa dS = OiS)+4T:{3Ro fi-—^+Oie)) = -47r+0(e)+0(<5) 

J\x\<3Ro J\x\=3Ro (3i?o)^ 

which is negative if e and S are chosen sufficiently small. This is a contradiction, 
and the claim follows. □ 
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5. Low ENERGY ESTIMATE 

Now we prove the low energy estimate in Proposition l3.1l which is the easiest of the 
three propositions to prove, especiaUy in the model case when there is no potential 
V, and the manifold M is a compact perturbation of R''. The idea here is to exploit 
the uncertainty principle to extract some gain from the spatial projection tp and the 
frequency projection Pjo, but in order to do this we need to somehow exploit the fact 
that H contains no resonances or bound states at the zero energy. For technical 
reasons (having to do with our use of the homogeneous norm iji/2(A/) instead 
of the inhomogeneous norm H^/'^{M)) we shall need the following non-standard 
formulation of this non-resonance property. 

Proposition 5.1 (Laplace-Bcltrami operators have no resonance). Suppose that 
f : M C is a measurable, weakly differentiable function such that 

f {xr''''^\f{x)\' + (x)-i-2-|v/|2 dg{x) < oo 

for some K < oo, and such that Hf = in the sense of distributions. Then (if 
a > is sufficiently small) f is a constant. 

Proof. We may take / to be real. The condition Hf = 0, combined with the 
local square-integrability of /, implies that / is in fact smooth thanks to elliptic 
regularity (and the smoothness of g). Since Hf ~ and H — Hq ~ outside 
of i?(0, i?,o), the function Aj^s/ is a smooth compactly supported function. Thus 
if we set F :— j^^^ * ^R^/i i-S- the convolution of Aj^a/ with the fundamental 
solution of the Euclidean Laplacian, then f — F is harmonic, F decays like 0(1/ (x)), 
and |VF| decays like 0(l/(x)^). From hypothesis and the triangle inequality we 
then have 

/ {xy'-^'lfix) - F{x)\' dg{x)< oo. 

JM 

Since f — F is harmonic, we conclude (e.g. from the mean-value theorem applied 
to f — F and its first derivative) that f — F is constant, thus f — F = c. By 
subtracting this constant from / we may in fact take F — f. This now shows that 
f{x) = 0{l/{x)) and Vf{x) = 0(1/ (x)'^). But then we can justify the computation 

/ \\/f\ldx^- f f{x)^^Vjf{x) dg{x)^ [ 2f{x)Hf{x)dg{x)=Q 

JM JM JM 

by inserting a suitable smooth cutoff to a large ball -6(0, R) and then letting R ^ oo; 
we omit the standard details. But then we have V/ = and we are done. □ 

Using this fact and a compactness argument, we now conclude 

Proposition 5.2 (Poincare-type inequality). Let f : M R be a Schwartz func- 
tion. Let ip be as in previous sections. Then for any e > we have 

(16) 

ll^|V/||U2(M) < Oe^o(l)(||(a;)-i/^--V/|U.(M) + ||(^)-''/'-^/||l^(m) 

+ £-i||(a;)-i/2--VH/|U.(M)+c-i||(x)-3/2-H/|U.(M)). 
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Proof. Suppose for contradiction that the Proposition was false. Then there exists 
a (5 > 0, a sequence e„ > converging to zero, and Schwartz functions /« such that 

+ e-i(:r)-i/2-Vi7/„|U.(M)+£-^||(:r)-3/2-i7/„||^.(^,)). 

Without loss of generality we may assume that /„ is not identically zero, and then 
we can normalize the expression in parentheses to equal 1. Thus 

||(:e)-i/2"-V/„|U2(m),||(x)-3/2-/„||^,(^^) < 1; 
(17) \\{xr'/^-^VHfjLHM) + \\{x)-^^^-''HU\\l2^m) < Sn, 

M^f\\\mM)>5. 

Next, we establish weighted bounds on /„ via a Bochner identity. From p7|l 
and Cauchy-Schwarz we have 

/ (a;)-i-2.Re((VVn)V~ffA:) < e„ = 0(1). 

We substitute H = — iV'^Vfe. Since M is flat outside of the compact set 5(0, i?o) 
we have 



V,F/„ = -iVjV'=Vfc/„ = -iv^V,Vfc/„ + 0(|V/„|)xs 
Using this and (|17|) . we obtain 

- / (a;}i-2.Re((VVrOV'=VjVfe^) < C. 

JM 



IM 

Integrating by parts we obtain 

/ (x)i-2-Re((v'=Wn)VjVfc7;) + (V'=(x)i-2-)Re((VVrOV'=V,Vfe7;) < C 

JM 

Since Re((VVn)^'^^j^*;/n) — ^Vfc|V/„|g, we can integrate by parts once more 
to obtain 

M ^ JM 

Since Am(2;)^^^'^ = 0((x)^^^^'^), the integral on the right-hand side is 0(1) by 
(ini). Thus 

(18) / (x)i-2-|Hess(/„)|2 <C. 

JM 

From this, H17|) and Rellich compactness we see that the sequence /„, when localized 
smoothly to any ball B{0,R), is contained in a compact subset of H^{B{0, Rj). 
From this and the usual Arzela-Ascoli diagonalization argument, we may extract a 
subsequence of /„ which converges locally in to some limit /. From (|17|l and 
Fatou's lemma we then see that 

(19) ||(x)-l/2--V/|U2(M), \\{x)-'/^-y\\mM) < 1; ||^|V/||U2(M) > s 

and that Hf = in the sense of distributions. But then by Proposition 15. II / is 
constant. But this contradicts the last estimate in H19|l . and we are done. □ 
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We can now quickly prove Proposition [3 .11 Applying lyTjl with uq replaced by PiqUq 
and HPioUo, and using some spectral theory to estimate the right-hand side, we 
obtain 

JO 

and 

JO 

note that Pio and H are spectral multipliers and hence commute with each other 
and with e~'*^. Applying Proposition 15. 21 Proposition l3.il follows. 



6. High energy estimate 



We now prove Proposition 13. HI which is the next easiest of the three propositions. 
This case resembles the local-in-time theory of Craig-Kappeler- Strauss [2| and Doi 
[S] , and indeed our main tool here will be the positive commutator method applied 
to a certain pseudo-differential operator, exploiting the non-trapping hypothesis to 
ensure that the symbol of the pseudo-differential operator increases along geodesic 
flow. As we shall now be working in the high energy setting, we will not need to take 
as much care with lower order terms as in previous sections. For similar reasons, 
we will not need to use the homogeneous Sobolev spaces H^{M), relying instead on 
the more standard (and more stable) inhomogeneous Sobolev spaces H''{M). The 
argument here is in fact quite general and would work on any asymptotically conic 
manifold with a short-range metric perturbation and a short-range potential. 

It will be convenient to use the scattering pseudo- differential calculus, which is an 
extension of the standard pseudo-differential calculus which keeps track of the decay 
of the symbol at infinity. We briefly summarize the relevant features of this calculus 
here, referring the reader to '5' for more complete details. For any m, ^ S R, we 
define a symbol a : T* M — > C o/ order (m, I) to be any smooth function obeying 
the bounds 

the function a(x,^) = {Cl"^ is a typical example of such a symbol. Note that 
we assume that each derivative in x gains a power of {x) , in contrast to the standard 
symbol calculus in which no such gain is assumed. We let S'™^'(M) denote the space 
of such symbols. Given any such symbol a G S''"''(M), we can define an associated 
pseudo-differential operator A = Op(a) by the usual Kohn-Nirenberg quantization 
formula 

Op(a)u(a;) {2ny" J e'^''-y^^^ a{x,^)u{y) dy d£.. 

We sometimes denote a by a (A) and refer to it as the symbol of A. Heuristically 
speaking, we have A = a{A){x, iV^). We refer to the class of pseudo-differential 
operators of order (m, /) as Also, if /i : R ^ C is any spectral symbol of order 

m/2, the corresponding spectral multiplier h{H) is a pseudo-differential operator 
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of order (m,0). In particular, (1 + i/)™/^ has order (m,0), and the Littlewood- 
Paley type operators Pio, Pmed, Pm have order (0, 0). We caution however that the 
Schrodinger propagators e~**^ are not pseudo-differential operators. 

The composition of an operator A = Op(a) of order (to, I) with an operator of 
B = Op(&) order (m', I') is an operator AB of order (m + to', I + /'), whose symbol 
a{AB) is equal to a{A)a{B) plus an error of order (to + to' — 1, Z + + 1); note the 
additional gain of 1 in the decay index I, which is not present in the classical calculus. 
Similarly, the commutator i[A, B] will be an operator of order (to + to' — 1,1 + 1' —1) 
with symbol <T{i[A, B]) equal to the Poisson bracket 

{aiA),a{B)} := VM^) ■ Vea(B) - V^^) ■ V,a(A), 

plus an error of order {m + m' — 2,1 + I' + 2). We shall write the above facts 
schematically as 

a{AB) = (t(A)ct(B)+0(5'"+"'-1''+''+1); a{i[A, B]) = {a{A), aiB)}+0{S'^+"'' -^'^+^'+^) 
or equivalently as 

Op(a) Op(6) = Op(a5)+0(*™+-'-i''+''+i), i[Op(a), Op(6)] = Op({a, 6})+0(*^+'"'-2''- 

In particular, since H has order (2,0) and has principal symbol ^1^1^(2;) plus lower 
order terms of order (1, 1) and (0, 2), we see that if a G S"^'\ then we have 

z[H,0p{a)]=0p{Xa) + 0i%!^/+'), 

where Xa denotes the derivative of a along geodesic flow in the cotagent bundle 
T*M. 

Associated with the scattering calculus are the weighted Sobolev spaces _ff™''(Af) 
defined (for instance) by 

Mnr^.^m \\{xy{l + H)"^/MLHM) 

(many other equivalent expressions for this norm exist, of course); when I = this 
corresponds to the usual Sobolev space H"'^{M). It is easy to verify that a scattering 
pseudo-differential operator of order (to, I) maps H™ (M) to iJ™ ~™'' +'(Af) for 
any m' ,1' . 

In (see also 0) it was shown that the non-trapping hypothesis on M allows 
one to construct a real-valued symbol a G 5^'" (depending on if) which was non- 
decreasing along geodesic flow, Xa > 0, and in fact obeyed the more quantitative 
estimate 

Xa{x,0^^{xml + \b\' 
for some symbol b of order (1, 1/2 — a). In Euclidean space, an example of such 
an symbol is C'v'{fy ' ^ some sufficiently large constant C'^. Quantizing this, we 
obtain 

i[H,A] = \7\{x)Vj + B*B + 0(*i^2-2<T-)^ 

where A := Op(a) is a symbol of order (1,0), and B := Op(&) is a symbol of order 
(1, 1/2 — a). We then apply the self-adjoint projection Phi — P^i to both sides, and 
observe that this commutes with H , to obtain 

i[H, P*hAPM\ = PmV^v{x)V,Pm + PLB*BPm + PMO{^lt^'')PM. 
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Applying (|14|l . and integrating by parts (discarding the positive term B*B, and 
using that '^l^''^" maps i/i/2 -i-h<t ^-i/2,i-<t^ obtain 

/ " C'II^/:»^^'»^IIhi/2(m)^h-i/^(a/)II"o|||i/2(m)+^^ ||Pwe-'*^Mo||^V2.-i+„ di. 

Since A is of order (1,0), and Phi maps the homogeneous Sobolev spaces to their 
inhomogeneous counterparts, we see that ll-F^i^^'/ij |liji/2(M)^^-i/2(A/) is bounded 
by some C^. To finish the proof of Proposition it thus suffices to show that 

||P/.,e~'*^^^o||^i/2,-i+. dt < o,,{l)\H0)\\h/2^My 
On the other hand, applying ^ to H^^ Phiu{Q) for j = 0, 1 we have 

^ ||Vi/-Je-"«P,„?/o||^o,-i/2-. + ||i/--'e-'*^Pw^^o|lHO,-3/2-„ dt < CK{T)^\\u{0)\\l,,,^^,j^ 

for j = 0, 1. Thus fixing t and setting / :— H^^e^^*^ PhiUo it will suffice to prove 
the fixed-time estimate 

1 

||-ff/||^i/2,-i+^(jVf) < Oeo^o(l)^ (^||Vi/-'/||^o,-i/2-<,(j\^) -I- ||i?VllH0,-3/2-„(M)^ 

It suffices to verify this for real- valued /. By Rellich compactness, the space 
iJi'-i/2-'^(M) embeds compactly into i/i/2'-i+'^(M). Since P^f ^ as Eq ^ 
for each individual /, we thus see by compactness that it suffices to show that 

ll'H"/llffl'-l/2-<7(M) < C'(l|Vi?/||^fO,-l/2-^(^,f) + \\Hf\\fjO,-3/2-,I^M) 
(20) + ||V/||^0.-l/2-^(JV/) + ||/||//0,-3/2-^(JV/)). 

The top order term of ||-ff/||Hi--i/2-<T(M) is already controlled by the right-hand 
side of l|2n|), so it suffices to control the lower order term ||i//||^o,-i/2-,T(^,/). But 
an integration by parts allows us to write 

Jm ^ Jm ^ ' 

The claim then follows from the Cauchy-Schwarz inequality. This concludes the 
proof of Proposition 



7. Medium energy estimate 



We now turn to the medium energy estimate, ProDOsition l3.2l which is the hardest 
of the three propositions. Neither the uncertainty principle nor the non-trapping 
condition will be much use here. Instead our tools^ will be a RAGE-type theorem, 
exploiting the fact that H has no embedded eigenvalues, to propagate the solution 
away from the origin after a long time r (though the energy localization shows the 
solution will not move too far away from the origin in bounded time), combined 
with a decomposition of phase space into incoming and outgoing waves (cf. Enss' 
method, [Hj, CHI); and several applications of Duhamel's formula. Because we 



One can also proceed via Kato's theory of _ff-smooth operators, using the Umiting absorption 
principles obtained in We will pursue this approach in detail in |2]. 
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have eliminated the high frequencies, we wiU enjoy an approximate finite speed 
of propagation law for the solution (but the upper bound of the speed is quite 

— 1/2 

large, being roughly 0(eQ )), and because we have eliminated the low energies, 
we will not encounter a frequency singularity when we decompose into incoming 
and outgoing waves, although again we will pick up some negative powers of Eq. We 
shall compensate for these Eq losses by using the RAGE theorem to gain a factor 
of 

Ot—>oo;€o,r{^) within a distance R from the origin, and to also gain a factor of 
o_R^oo;eo(l) when one is further than R from the origin. The reader should view the 
comparative magnitudes of eo,R,T according to the relationship I/eq ^ i? ^ t, 
which is of course the most interesting case of Proposition 13. 21 

It is instructive at this point to recall the basic features of the Enss' method. Define 
the wave operator W+ = s — limt_>+oo e^*^e~^*'^° . Enss' method is concerned with 
establishing completeness of VF+, i.e., showing that the range of W+ coincides with 
the continuous subspace oi H = Hq + V. We assume otherwise so that there exists 
(po not in the range of W+ . Density arguments allow us to consider (/jq with compact 
support and medium energies only. We then evolve 4>o by e~**^ and claim that we 
can find a sequence of times tn and a decomposition 

(21) e''*"-^(/;>o = 4'n = 4>n,out + 0n,m + 4>n, 

(compare this with the decomposition into Fioc, Fgiob in H24(l and a further decom- 
position of Fgiob into the outgoing/incoming waves in (PU) with the properties 
that 

||(W+ - l)4>n^out\\L^, ||0n,m||L2, ll^nlU^ = 0„_>oo(l)- 

If we had such a decomposition, we could conclude from the boundedness of W+ 
that 

and thus (po lies in the range of W+ . The desired decomposition H21|l can be found 
for instance in JHIi and we sketch it as follows. The local component </)„ can be 
set for instance to (/>„ := X\x\<n'Pn, its convergence to zero is a consequence of the 
RAGE theorem (if the times i„ are chosen appropriately). The global component 
0'~X\x\<n)4'n is partitioned into functions (f)n,a supported in unit balls with centers 
at the lattice points a G R". One can then define the 

^n,a,out^ e^^''~y^-^m[Ci'l>nAy)dS.dv, (t>n,a,zn^ e^^''^^^'^ ( l-m(0 )</>„,« (?/) c?2/, 

where m{^) is a smooth multiplier localizing to the region Z(^, a) < 37r/2. Since 
it is expected that a compact support function 0o propagated by e~'*"-^ should 
become mostly outgoing in the region > n we have that 0, while the 

fact that (W+ — l)(j)n,out — > in follows since on the outgoing waves the evolution 
g-ttH .^gj]^ approximated by the free flow e~'^*^° and in fact converges to it as 
n — !■ oo. 

Our proof of Proposition 13 . 21 will follow in spirit the above construction, although 
we have to take additional care since we do not work with functions of compact 
support and need weighted estimates instead of bounds. On the other hand, 
we are still localized in energy, and so we will not be too concerned about losing 
or gaining too many derivatives (as we are able to lose factors of Eq or Eq^ in our 
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estimates here). In particular, the metric perturbation H of Hq is now of similar 
"strength" to a potential perturbation, and we will now be able to use Hq as a 
reasonable approximant to H, at least when the solution is far away from the 
origin. 

It is more convenient to work in the dual formulation. First observe that the claim 
is easy when T < r, since Pmed^~^*^ maps Fi/2(M) to H^{M) thanks to the 
frequency localization of Pmed- Similarly the /J" portion of the integral is easy to 
deal with. Thus we may assume that T > r and reduce to proving 

^ ||¥'VP„ede-'*^no||i2(M) dt < (Ceo,iJ,r+Ofl^cx>;eo(^(2^))+Or^cx>;eo,ii(i^(2^)))ll^*o|||i/2(M)' 

which after dualization (and shifting t by r) becomes 
(22) 

(Kim 

{f ^ \\F{t + T)\\l.^^^dtY'^ 
Jo 

for any vector field F{t) defined on the time interval [r, T] which is Schwartz for 
each time t. 

Fix F. We split 

(23) e''"e'-"PmedVj^F^{t + r) = e''"Fioc{t) + Fg^it) 

where 
(24) 

Fioc{t) := ^Re'^"P^edyj^F^{t + T); Fgiob{t) := {1 - <pR)e'^" P^edVj'pF^t + r). 
• Term 1: Contribution of the local part. 

Consider first the contribution of the local term Fiodt). To control this term we 
use the following local decay result. 

Proposition 7.1 (RAGE theorem). Let (pn he a biirnp functAon supported on 
B(0,2R) which equals 1 on B{0,R). For all Schwartz vector fields f^, we have 

Pmed'^ j'Pf \\l^(M) < OT^oo;R,eo{\\f\\L^(M))- 

Proof. Observe that we have the crude bound 

\\ipRe'''"Pmed'^j'PfmLHM) < CsJf\\H-i.-^M) 

since Pmed^^ maps H~^'~^ to L^, and (Prb^'^^ is bounded on L^. Since the unit 
ball of L^{M) is precompact in H~^'~^, it thus suffices to prove the estimate 

\\(PRe'''^" Pmed'^ j^Pf h'^iM) < Or^oo;fl,/,£o(l) 

for all Schwartz vector fields /; the point being that the compactness allows us to 
ignore the dependence of the constants on /. Fix /, (f, R. Since / is Schwartz, the 
curve {'pe^'^^ Pmed^ j^f'' '■ T G R} is bounded in H^'^{M) (for instance) and hence 
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precompact in . Thus to prove the above strong convergence, it wiU suffice to 
prove the weak convergence result 

Or^oo;fi,/,eo,i/.(l) 

for all Schwartz functions ■0. 

Fix ■)/). Since the spectrum of H is purely absolutely continuous'^, we see from the 
Riemann-Lebesgue lemma and the spectral theorem that 

and the claim follows. □ 



From this proposition we sec in particular that 

\\Floc{t)\\LHM) < Or^oo;fl,£o(l)ll^(* + T)llL2(M)- 

From dualizing the second part of ITJ we have 



(25) II / e^*^(x)-^/^--G(t) di||^_,/.(^,) < K{T){ / ||G(i)||i.(M) dt) 

<J <J 

for any G. If we truncate the time interval to T— t, substitute G{t) :— {x}'^/^'^'^ Fioc{t), 
and take advantage of the spatial localization of Fioc to 5(0, 2R), we thus have 

II / e''"Fiocit) dt||^-i/.(M) < o.^oo;e„,fl(if (r))( / ||F(i + T)||i.(^,) 
JO Jo 

Thus the first term in (|23|l is acceptable. 

• Term 2: Contribution of the global part. 

To conclude the proof of Proposition 13. 21 it suffices to establish the estimate 
(26) 

e''"Fgiob{t) dt||^-V2(Af) < {Ce„,R,r+OB^^oo;eoiK{Tm / ||F(t+r) || 



for the global term Fgiob{t). 

Let us first show an associated estimate for the free propagator e'*^" , namely 

(27) II r \''"'>Fgiob{t) rft||^-v2(M) < Ce„,RA f ^ mt + r)\\l.(M) dt)'/'. 
Jo JO 

Note that our constants are allowed to depend on r, as we will no longer need to 

place a factor of K{T) on the right-hand side. 



To prove (|27|l . first observe that by dualizing the second part of we have 



e 



■^Actually, the argument would still work well if H had some singular continuous spectrum, 
except that r must now be averaged over an interval such as [tq, 2ro], but the reader may verify 
that the arguments below will continue to work with this averaging. The spectral fact which is 
really being used here is that H contains no embedded eigenfunctions at medium frequencies, 
since such eigenfunctions would certainly contradict i6l . 
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(for instance), so it will suffice to show that 

\\Fglob{t)\\H°'^(M) < Cso,B.,T\\F{t + t)\\l2(^j^,j-^ 

for all < t < T — r. On the other hand, we know from inspection of the symbol 
of Pme(Nj<P that 

||P™edV,</.F^(t + T)||H2,2(M) < C,jF{t + T)\\L-2(M) 

(for instance), so it will suffice by 124|l to show that 

But this can be easily established by standard energy methods''. This proves (|27|) . 
Thus to prove H26|l it suffices to show that 
(28)^ , 

II ( '(e^*^-e^*^o)F,i<,fc(t)dt||^„,/.(,,,) <Ofl^oo;eo(l+if(r))(/ '||i^(t+r)||i.(M) 



At this stage it is necessary to decompose Fgiob{t) further into "incoming" and 
"outgoing" components, which roughly correspond to the regions of phase space 
where x • ^ < and a; • ^ > respectively. Semiclassically, we expect Fgiob{t) to 
be supported almost entirely in the "incoming" region of phase space, since it is 
currently far away from the origin, but came by propagating a localized function 
backwards in time from t+r. However, if this function is supported in the incoming 
region of phase space, then by moving further backwards in time by t it should move 
even further away from the origin, and in particular it should evolve much like the 
Euclidean flow (i.e. it should become small when e**^ — e**-^" is applied). 

We now make this intuition precise. The first step is to formalize the decomposition 
into incoming and outgoing waves. We first take advantage of the fact that the 
spectral support of Pmed vanishes near zero. From ^ we have 

(29) Fgiobit) - (1 - ^B.)H^Fgiob{t) 
where 

(30) Fgiobit) e'^"P,-aedy,^F\t + r) 

and Pmed ■= H~'^Pmed- This factor of we have extracted from Pmed shall 
be helpful for managing the very low frequencies in the proof of Proposition 17.21 
below, which would otherwise cause a significant problem for this portion of the 
evolution (at least in three dimensions; this step appears to be unnecessary in five 
and higher dimensions, for reasons similar to why resonances do not occur in those 
dimensions). 



We now need the following phase space decomposition associated to the Euclidean 
fiow e-**^o. 



"'For instance, if / g H'^'^{M), one can first estabUsh uniform bounds on e^'^^ f in H'^'^(M) 
by spectral methods, then use energy methods to control e'^^/ in H^'^(M), and then finally 
in H^'^{M), losing polynomial factors of r in each case. One could also argue using positive 
commutator methods based on I13i for such operators as A= (x)*. We omit the details. 
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Proposition 7.2 (Phase space decomposition). There exist operators Pin,Pout 
such that 

(31) (1 - ipR)H' = (1 - ipR)H'P,„ + (1 - ipR)H^Pout 
and for which we have the estimates^ 

(32) ||e-^o(l - ^R)H^P,„f\\H2.-s^M) < CiR" + |s|)-i-- ||/||^.o(,,) 
and 

(33) ||(x)3/2+-(l _ cpR)H^P,^te^^"«fU.^M) < C{R^ + \s\)-^-^\\f\\m^.^^(M) 
for any time s > and all Schwartz f . 

The proof of this proposition is a straightforward application of the principle of 
stationary phase, and can be justified heuristically by appealing to the intuition of 
microlocal analysis and the uncertainty principle. It is however, a little technical 
and will be deferred to the next section. Assuming it for the moment, let us conclude 
the proof of Proposition 13. 21 It suffices to prove (|26|l . From H29|l and Proposition 
13. 21 we can split Fgiob as 

(34) Fgiob{t) - (1 - VR)H^P^nFglob{t) + (1 - ^ r)H^ PoutFglob{t) , 

and treat the components separately. 

• Term 2(a): Contribution of the global incoming part 

We now control the contribution of the incoming component of to (jJEJ- We 
use Duhamel's formula to write 

{e''"-e''""){l^^R)H^P,nFgiob(t) = z f e'^'-'^" {H-H^)e''"<> (1-^r)H^ P.r^Fgiobit) 

Jo 



and so it suffices to show that 

\£ ^ j\'^*~'^"{H-Ho)e''""{l-VR)H^P,nFgiob{t) dsdiW^-i,,^^^ 



= OR^^,,,{K{m / \\F{t + T)\\l,,^,. dtf'^ 







Substituting t' :— t — s using Minkowski's inequality, we can estimate 

1/ / e<'-''^"iH-Ho)e''"»il-^R)H^P,„Fgiobit) dsdt\\^^„,,^,j. 
Jo Jo 

< r ^\\ r ^ ' e'''"{H - Ho)e''"°P,nFgiob{t' + s) dt'\\^_,,,,j^,jy 
Jo Jo 



^The decay weights here should not be taken too seriously: indeed since we are assuming H 
to be a compactly supported perturbation of Ho we have enormous flexibility with these weights. 
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Applying if^ . it thus suffices to show that 
(35) 

I IP/4- 1 _M|2 

On the other hand, since H — Hq is a compactly supported second order operator 
we have 

\\{xf^^+^H-Ho)e''"'{l-VB)H^P^nFglob{t' + s)\\L^M)<C^^^^^ 

and so from H32|l we thus have 

\\{xf^'+''{H-Ho)e'''"°{l-VR)H'P.nFgUt'+s)\\L2^M)<C{R'+\s\^^^^ 

From 1(201) we note that \\Fgiob{t' + s)\\l2(m) < ^£011^111,2(4/) (noting that P,ned 
maps H^^{M) to H'^^{M) with an operator norm of C^^. Thus we have 

r " '||(x)3/2+-(i7_i7o)e-ffo(i„^^)^2p^^^^^^^^(^,^^)||^^^^^^ 
Jo 

l-T-T 

i-'^r / II pr^ 4- ^^11? ^^^1/2 



<C,,{R' + \s\r'~% lini + T)||i.(M) dt) 

Jo 

and the claim H35|l follows upon integrating in s. 

• Term 2(b): Contribution of the global outgoing part 

We now control the contribution of the outgoing component of (|;-{4|l to (|26|l . From 
(|25|l we have 

II r ^ e''"{l~^n)H'P,^tFgiob{t) dt|lH-i/^(M) < K{T){ C ^ \\S{t)\\l.,^,. dtf'^ 



2{M) ^ ^M-t n / 11-^ V'-y II L2(M) 

j\) JO 

where 

S{t) := (x)3/2+-(l - ^R)H^PoutFglob{t). 

From (0J, a similar estimate holds when e^^^ is replaced by the free flow e'*^" , with 
K{T) replaced by a constant C. Thus to control this contribution to (|26|l it will 
suffice to show that 

li^(M) dtf'^ = o«^o;.o( f ^ ll^^(i + r)||i.(M) dtY/\ 

JO 

To prove this, we first use (|3()|l to expand 

l|5(i)lli^(A/) - ((a;)'/'+"(l - (^i?ji/2Ponte'"''PmedV,(pF^(i + t),5(<)) 

= -{F\t + t), vpV,P,„ede-'^^T4^(t)) 

<C||F(i + r)|U2(M)||P,„ede-'"''M/(t)||ffi,-20(M) 

<C,jF{t + T)\\L2i^M)\\e-""W{t)\\H-^o,-.o^M) 

by the support of Lp and the smoothing properties of Pmed, where 
(38) W{t) P:^,H^{1 - ^«)(x)=^/2+-5(t). 
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We use Duhamel's formula to write 

e~'^"W{t) = e-'^"°W{t) - i f e-'^'-'^^iH - Ho)e-''"'>W{t) ds, 

Jo 

and apply iJ~20,-20 ^Qj-j-^g on both sides, to obtain 

\\e-^^"Wit)\\H-2o.-.o^M) < \\e-^^"^'Wit)\\H-2o.-2o^M)+ r \\e-^^^''^"{H-Ho)e~^'''«W{t) ds\\H-.o.-.o^My 



Observe that the propagator e^*'^'^^'*'^ maps H^^^{M) to H^^^{M) and hence to 
^-20,-20^^^^^ uniformly in t and s. Also, since H — Hq is a compactly supported 
operator we see that H — Hq maps H^^^'~'^^{M) to H^^"{AI). We thus have 

\\e-""wmH^Bio,c)) < C\\e-""«W{t)\\H-^s.-..^,j^+C f || ||e-^^^°M^(t)||^-i8,-ii(M) ds, 

Jo 

and thus by (|38|l and the adjoint of (|33|) 

\\e-'^"W{t)\\HiiBio,c)) <C{R^ + T)-^-^\\S{t)\\L2^M) + C r {R^ + s)-^--\\S{t)\\L2(M) 

Jo 

<CeoR^''\\S{t)\\L2(^M)- 

The point here is that we have obtained a non-trivial decay in R. Inserting this 
estimate back into ()37(l we see that 

\\SmLHM)<CeoR-'\\F{t + r)\\L2(M), 

and H36|l follows. 

The proof of Theorem II. II is now complete, once we complete the proof of Propo- 
sition which we do in the next section. 

8. Proof of the phase space decomposition 

We now prove Proposition 17. 21 The ideas here have some similarity with a decom- 
position used in ^2 ; related ideas were also used recently in [31 • The idea of using 
a decomposition into incoming and outgoing waves to analyze the perturbation 
theory of the free Laplacian Hq of course goes back to Enss (and, in a different 
context, even earlier to the work of Lax-Phillips.) The phase space decomposition, 
developed by Enss |S] and refined by Simon [TO], is based on the following con- 
struction. Let {Cj}jg2'' ^^"^ ^"^jJ^jeZ'' smooth partitions of unity in and 
R| respectively with the property that each of the functions Cj{x) and mj{£_) is 
supported in the ball B{i, 2). Define the symbols (of the \1/D0's) of the projections 
Pin and Pout on the incoming and outgoing states: 

Ptn = ^ mj(^)Ck(a;), Pout = ^ TOj(OCk(a;) 

jk<0 jk<0 

The following estimate, crucial to the Enss' method, reflects the expectation that 
the outgoing waves never come back to the region where they originate. For any 
t > 0, N > and j • k > such that |j | > 3, we have 

l|e-*"''™jCk/||,.(,^,<,/,,„,,/,) < C.d + t + |k|)--||/||,,(R.^ 
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(compare with H32() noting that the transformation t —t eorresponds to the 
change Pin — > Pout)- Observe that the condition |j| > 3 ensures that we only deal 
with the outgoing waves of velocities bounded away from zero in absolute value. 

A continuous decomposition, based on coherent (Gaussian) states, was used by 
Davies, while the outgoing/incoming waves defined via projections on posi- 
tive/negative spectral subspace of the dilation operator x-V + V -x were introduced 
by Mourre, 0. 

We now return to our decomposition into the incoming/outgoing wave needed to 
prove ProDOsition l7.2l We first observe from dyadic decomposition that it suffices to 
prove the claim with (1 — ipfj) replaced by {ip2R — Vi?), since the original claim then 
follows by replacing R by 2"^R and summing the telescoping series over to > 0. By 
further decomposition, and some rotation and scaling, we may replace f2R—VR by a 
smooth cutoff function -0 = ifjn supported on the ball B := {x : \x — Rei\ < i?/100}, 
where ei is a basis vector of R". 

Our construction will be based on Euclidean tools such as the Fourier transform, 
so we shall now use the Euclidean inner product and Euclidean Lebesgue measure 
instead of the counterparts corresponding to the metric g. In particular the operator 
H will no longer be self-adjoint, but this will not concern us as we shall soon break 
it down into components anyway. To reflect this change of perspective we shall 
write our manifold M now as R". 

We begin with the Fourier inversion formula 

Jr" JR" 

valid for all Schwartz /, where -0 is a smooth cutoff to the ball B := {x : — i?ei| < 
R/50} which equals 1 on B. We then split the ^ integration into subspaces > 0, 
defining 

J{i<o JR" 

and 

Poutfix) = / / e2-(--«H^(y)/(zj) dyd^ 
Jii>a JR." 

where := ^ • ei is the ei component of ^i. We remark that these operators 
are essentially Hilbert transforms in the ei direction; the multiplier is of course 
discontinuous in the variable, but we will never integrate by parts in this variable 
so this will not be a difficulty. 

Clearly we have the decomposition (|^ . It remains to prove 
• Proof of (ESI 



We now prove the estimate ()33|l . Since H = Hq on the support of ?/;, we may write 
{xf/^+'^ijH^ = i?-^/2+"a(x)(i?V,)4 
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for some bounded functions tensor a{x) supported on B. We shall think of RVx as 
a normalized gradient on B. 

In light of the above decomposition, it thus suffices to show that 

J >o ^R. 
for all all times s > 0. 

We first dispose of the derivatives (i?Vj:)*. Each x derivative in i?V^ hits the phase 
^2-Ki{x-y)-i ^ where it can be converted to a y derivative, which after integration by 
parts either hits or e^''^°f{y). Since partial derivatives commute with e**^", we 
thus reduce to showing that 
(39) 

f e'^^^^-yy^i.^^^Hy)e-"°f^^^Hy)dydaLUB)<CR^~^^ 

whenever < ji,j2 with ji + j2 = 4, where -i/)'-'^' := {RVxy^-fp is a minor variant 
of the cutoff 1^, and /(j'^) := {RV f ■ 

Applying a smooth cutoff to /, we can divide into two cases: the local case where 
/ is supported on the ball {|a;| < i?/100}, or the global case where / is supported 
on the exterior ball {x > i?/200}. Let us consider the global case first, which is 
rather easy and for which one can be somewhat careless with powers of R. We first 
observe from Plancherel (or the boundedness of the Hilbert transform) that 

Thus to prove H39|l in the global case it suffices to show that 

We divide further into two sub-cases, the short-time case s < R^ and the long-time 
case s > R^. In the short-time case the claim follows since 

= ^II/'"^IL^(R^) 

which is certainly acceptable (if a is small enough). In the long time case we 
interpolate between the two bounds 

ll^^"^^'^''"5lL.(R3) <C||e-^05lL.(R3) <C||5l|,.(R3^ 

and 

W^^'^^e-^OgW^^^^s^ < Ci?"/2||e-^°g||^^(R3^ < CR'/'s-'/^g\\^,^j^;.^ < Ci?"/2.-"/2||5||^o,„,.,„(R3^ 
to obtain 

||^Oi)e-i/ojfc)||^^^^3^ < C(i?/,s)-i-'^||/(^=)|iHO.+c. < C(i?/s)-i-'^i?^"^i?-"+i+^ni/ll^.s,.i(R3) 
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which is certainly acceptable. 

It remains to prove (|39|l in the local case, when / is supported on the ball {|a;| < 
i?/100}. In this case expand the fundamental solution of e^^^° to write 

where 

Observe from all the spatial cutoffs that x, y are localized to the ball B, while z is 
localized to the ball \z\ < i?/100. Also, ^ is localized to the half-plane > 0. Our 
task is to show that 

(40) II / K,ix,z)f(^^Hz) dzU^B) < Ci?^-^(i?' + s)-'-ni/ll^is,ii(R^- 

J\x\<R/100 ^ ' 

We split Kg into and Kl°, corresponding to the regions |^| > and 
1^1 < R~^~^'^ of frequency space respectively. The contribution of K!}'^ will be very 
small. Indeed, for any |f | > _R^^+°', we can evaluate the y integral using stationary 
phase as follows. Observe that the y gradient of the phase 

2tt{x - y) ■ ^ - \y - z\y2s 

is equal to 

-27rf - (y - z)/s. 

From the localizations on y, z, and ^ we observe that this quantity has magnitude 
at least > c(|^| + R/s). One can then do repeated integration by parts in the yi 
variable, gaining an R every time one differentiates the ip^^^^ cutoff, to obtain the 
bound 

I f g2^»(x-y)-?^^o-i)(y)g-»|a-2|V2. dy\ < CNR'^iRl^l+RVs)-^ 

for any N. Integrating over all |^| > we thus sec that K'^'^{x, z) is bounded 

by Om{{R^ + s)^^) for any iV, and so this contribution to (|40|) is easily shown to 
be acceptable (using crude estimates on Z'--'^-'). 

Now we deal with the low frequency case |^| < R^^^" . In this case we expand f'^^^^ 
and integrate by parts to write 

K[°{x,z)f^^^Hz) dz\ < CR^'l [ Vi'Ki°{x,z)f{z) dz\ 

z\<R/100 J\z\<R/100 

< Cs-^'^ , 

' \z\<R/WQ J \i\<R-'^+' ^iiyQ 

e2''*(^-^)-«V'^''Hy)e"''^"^''/^' dy\W^\^\i''\f{z)\ d^dz. 

/R' 

We crudely bound i^-'^l^p^ by 0{R'^'^). From stationary phase we also observe the 
estimate 

I / e^M^-y)<^Ui)(^yy-^y-^\V2s < ^^3/2 

7r' 
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while by taking absolute values everywhere we also have the estimate 
and hence 

K[°{x,z)f^'^\z) dz\<Cs-'^/^R^" f I min(s"/2,i?3)|/(z)| d^dz 

2|<i?yioo "'|z|<fl/ioo J|ci<i?-i+<' 

< Ci?'^'"s-3/2j^in(s3/2^i?3)^-3 j 

J\z\<R/WO 

<Ci?^^s-3/2min(.3/2,i?3)i?-3/2||/||^,(j^3^ 

which is certainly acceptable. This completes the proof of H33|l. 
• Proof of (p!^ . 

To conclude the proof of Proposition 17.21 it suffices to prove H32|) (with 1 ~ ipj^_ 
replaced by "0, of course). This is in a similar spirit to the proof of H33() . although 
the steps will be in a somewhat permuted order. We begin by estimating 



fe=0 
2 

k=0 

Now observe that 

4 

for some smooth cutoff functions aj^k adapted to the ball B. Our task is thus to 
show that 

||(x)-«e-^°a,, feVi+^-p„/||^,(j^3^ < CR'-^{R' + ||/||^,,(R3^ 
for aU < j < 4 and < fc < 2. 

Fix J, k. Let us first control the contribution of the weight arising from the 

region \x\ > R/100. This is contribution dealt with differently in the short-time 
case s < R? and the long-time case s > R^. In the short-time case s < R^ we 
control this contribution can be controlled by 

<Ci?-8||e-^0a^-fcVJ"+'=P„/||^,(j^3^. 

Since e*^^" and aj^k are both bounded on L^, we can bound this by 

<Ci?-«||VJ"+^P™/||^,(P^3j. 

Now observe from Plancherel that Pin is bounded on H^'^ (since multiplication by 
tp is certainly bounded on and so we can bound this by O(i?^^||/||//2o), which 
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is acceptable. In the long-time case, we control the contribution instead by 



(41) < Ci?-'+^||e'^^V,,feV^+'=P™/||^„,_3/.(R3 



Now we interpolate between the energy estimate 
and the decay estimate 
to obtain 

The operator aj^k maps to H'^^^/'^ with a bound of 0(i?'^/^), so we can therefore 
bound (dH) by 

As in the long-term case we can bound ll^"'^''-Pjn/ll^2(R3j by ||/||^20(R3)j and so 
this case is also acceptable. 

It remains to control the contribution in the region |a::| < i?/100. We then expand 
out the fundamental solution of e*"^" , and reduce to showing that 

(42) 

Note that the a; variable is localized to the ball \x\ < i?/100, while y and z are 
localized to the ball B, and ^ is localized to the half-space < 0. Once again, we 
split into the high frequencies |^| > and low frequencies |^| < R~^^"'. In the 

case of the high frequencies, we move the y derivatives onto ip{z)f(z) by integration 
by parts, and reduce to showing that 

S-3/2||(a;)-8 /" f ^-i\x-y\'/2s (^yy■.^(y^z)■i 

7|«|>i?-l+<';6<0 JR' JR" 

g{z) dzdC|U2(|,|<fl,/ioo) < CR^-^{R^ + \s\)-^-"\\g\\^.^s) 
for some g. But the y derivative of the phase 

-|x-y|V2s + 2^(z/-z)-e 

is 

-{x-y)/s + 2TT£, 

which has magnitude at least c(|^| -I- R/s), by the localizations on x, j/,^. Thus by 
stationary phase we have 

e-'l^'^l'/2-aj,fe(y)e2"(^--)-« dy = ON{R'm + R/s)-^) 

R' 

for any > 0, and the claim is now easy to establish by crude estimates. Thus 
it remains to prove (|42|l in the low frequency case \^\ < In this case we 
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convert the V^+'' derivative to 0{\£,\^+'') = 0{R-^^+''^ R^"), and thus estimate the 
left-hand side of this contribution to (|42|l by 

JTi, 

But by stationary phase as before, we have the estimates 



sup 

|4|<fl-l+<';{l<0 ' JR^ 



/R 

and thus we can bound the previous expression by 

<Cs-y^R-^^+'^^R^^R-^mm{s'^',R'm^r' f .^V'WI/WI d2||L=(|.|<fl/ioo)- 

^R 

We crudely bound the L'^{\x\ < i?/100) norm of (x) ^ by 0(1), and use Cauchy- 
Schwarz we can bound the previous expression by 

< Cs-3/2i?-(j'+'=)i?C-i?-3min(//^i?3)^3/2||^||^^^^^^^ 

which is acceptable (treating the s < R^ and s > R^ cases separately). The proof 
of Proposition 17 . 21 is now complete. ■ 
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